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Abstract—An explicit, near optimal, closed-loop steering logic common to both, 2- and 3-dimensional 
launch trajectories is presented, This scheme uses a single steering angle in a canted plane, which can then 
be resolved as pitch and yaw angles in any inertial frame. New, exact formulations for integrating the 
steering thrust acceleration are presented along with the existing series approximation. Using Enke’s 
method, two uniform gravity vectors are evaluated to account for the actual effects of inverse square 
gravitational force field in the propagation of position and velocity vectors, during the future guided 
trajectory. An iterative differential corrector to solve for the guidance parameters that decide the steering 
angles and thrust cut-off is discussed. Also, means of overcoming the difficulties faced by the closed-loop 
guidance logic close to injection is presented. Simulations are conducted for a typical launch mission. The 
computational load for the onboard computer are estimated in order to explore the feasiblity of realizing 
this guidance logic. Simulation study of the effects of computational delay is also presented. 


I. INTRODUCTION 

Satellite launch vehicles are characterized by many 
uncertainties due to rapid burning of fuel, swift 
changes in vehicle parameters, high accelerations, 
discontinuous thrusting of multistage vehicles and 
changes in the environment. To use such vehicles for 
accurate orbital injection of the payloads, a highly 
accurate, optimal, closed-loop guidance (CLG) logic 
is required. Over the last three decades many guid¬ 
ance schemes have been developed [1-3]. The explicit 
methods determine in real time the steering com¬ 
mands, based on present and target states. Most of 
these schemes employ certain optimization criteria 
for improved performance. Since, real time solutions 
are essential for explicit guidance schemes, simplify¬ 
ing assumptions are necessary to get analytical sol¬ 
utions for the highly non-linear two point boundary 
value problem. At each update cycle, the parameters 
of optimal logic are calculated depending on the 
current and the target conditions only. The repeated 
updates of guidance parameters can ensure the attain¬ 
ment of the terminal conditions through redefinition 
of the optimal trajectory at each cycle. Therefore, the 
explicit schemes are flexible and capable of overcom¬ 
ing larger uncertainties. They permit redefining 
of target conditions during the flight itself for the 
mission salvage. They can also be easily adapted 
for different missions with the same rocket. Hence, 
the explicit guidance schemes are preferable for 
developmental vehicles where the uncertainties are 
significant. 

Developmental efforts in the explicit guidance 
schemes can be traced back to Lawden’s [4] deri¬ 
vation of optimal thrust direction in an analytical 
form for, uniform gravity, constant thrust and mass 


flow rate and vacuum conditions of flight. The opti¬ 
mal thrust vector orientation is given by the bilinear 
tangent law of the form: 

+Bt)/(C+Dt), 

where is the thrust attitude angle. A, B, C, D are 
the steering coefficients determined from the bound¬ 
ary conditions and t represents the time variable. 
When the horizontal range is a free terminal par¬ 
ameter, while terminal velocity and altitude are 
specified, the coefficient D in the above bilinear 
tangent law becomes zero, leading to the linear 
tangent law as: 

tanfi/'jt)] = A' + C't. 

Chandler and Smith [5] extend the linear tangent law 
to develop a 3-dimensional guidance scheme. The 
steering logic for the yaw and pitch are given in the 
decoupled form, under the assumption of uniform 
gravity. The uniform gravity assumption with the 
equivalent fiat Earth approximation is valid for a 
short launch trajectory with small magnitudes of 
vehicle velocity. 

Cherry [6] overcomes the limitations due to the 
uniform gravity field assumption, by judiciously com¬ 
bining the sphericity of earth with the guidance law, 
as: 

sin[+(0] = A + Bt - (g • y)/a T , 

sinW'p(t)] = C + Dt - (p/R 2 + Vji/R)/a T , 

where <l/ p (t) are the thrust attitude angles in 
yaw and pitch directions, p is the gravitational con¬ 
stant, a T is the thrust acceleration magnitude, V H is 
the vehicle velocity along horizontal direction in orbit 
plane, R is the radial distance of the vehicle projected 
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on to the orbit plane, y is the position vector perpen¬ 
dicular to the orbit plane, g is the acceleration due to 
gravity. Teren [7] uses similar method for removing 
gravity effects entering into the dynamical equations. 
It is observed that, as the target is approached, 
attitude turn rates toward instability due to the 
numerical errors and requirement of zero injection 
error to be met in the remaining short time. Teren [7] 
suggests termination of CLG calculations except for 
Tgo parameter (time to thrust cut off), near the final 
injection. The thrusting is terminated on attaining the 
required injection energy level. 

Bittner [8] derives a tangent form of steering law by 
taking into account the coupling between pitch and 
yaw as shown below: 


tan[i/r y (r)] = A + Bt, 
tan[t/t p (0] = (C + Z)r)cos[i/t y (/)]. 


He suggests the minimization of the performance 
index in the form of terminal error and total burn 
time instead of demanding exact injection as. 


y = (d-x T ne](d-x T )- 


d/, 


where d is the desired terminal state, X T is the 
estimated injection state and [Q] is a suitable weigh¬ 
ing matrix. Thus, the numerical instability close to 
injection is avoided. The inflight parameter identifi¬ 
cation is recommended to ease out instability due to 
the off-nominal performance of launch vehicle. After 
studying the effects of truncation in on-board com¬ 
puter, Bittner concludes that a 24 bits fixed point 
arithmetic is necessary. 

Lu [9] derives a linear tangent steering law, similar 
to Marec[10] as: 


4(0 = (X + ku - T ; Mi + ;?(t - 


where A f (f) is the instantaneous unit thrust direction, 
X is the unit reference thrust direction, T, is the thrust 
vector time constant. For 3-dimensional steering also, 
the vectors X, A define a canted plane containing the 
thrust vector. The pitch and yaw angles are obtained 
from single steering angle within the above canted 
plane. To get better accuracy, equivalent average 
gravity vectors Saw » 8avr are separately defined to 
represent the effects of gravity on estimates of vel¬ 
ocity and position vectors respectively. 


Saw = (1/A)| X g, - (1/2) (go + g s 

J-o 


Savr — (\!N 2 )\ 2 X (N -j)gj-(N + l/3)g 0 + (l/3)g. 

. j = o 


*] 


where g y is the gravity vector at the j th segment on the 
future trajectory. 

Recently, Sinha et al. [11-14] developed a very 
elegant explicit guidance scheme keeping in view 
requirements of the Polar Satellite launch vehicle 
(PSLV). The average gravity vectors are computed 


using the Enke’s method. The guidance scheme has 
been extended to consider non-uniform mass flow 
rate and thrusting which are typical of solid propel¬ 
lant engines. 

In this paper, further studies on this guidance 
scheme are attempted. Analytical closed-form evalu¬ 
ation of thrust effects are proposed against approxi¬ 
mate evaluations given earlier. Typical on-board 
computational time requirements are studied. 

2. OPTIMAL GUIDANCE LOGIC 

The objective of the optimal guidance logic is to 
determine the thrust attitude angle by a closed-loop 
action, such that the multistage launch vehicle places 
the payload/satellite into the desired orbit with mini¬ 
mum thrusting time (time-to-go). The guided trajec¬ 
tory of the vehicle is truly 3-dimensional, since, the 
plane containing vehicle trajectory at the launch 
point and that of the final orbit are non-coplanar. 
The algorithm is tested for the PSLV-class of vehicles. 

Launch Point Inertial (LPI) co-ordinates and Orbit 
Plane Geocentric Inertial (OGI) co-ordinates are used 
in the analysis. In the LPI co-ordinate system, the 
axis /, is along the radius joining earth center to the 
launch point. I x is parallel to launch azimuth and /„ 
completes the right-handed co-ordinate system. The 
OGI reference frame, shown in Fig. 1, is used to 
estimate a suitable injection point on the orbit. The 
five orbital parameters, a, e, i, (i and 9 are specified 
while leaving argument of perigee (to) at injection as 
an unspecified parameter. Then, is defined along 
the projection of the present position vector of the 
launch vehicle, on the orbit plane and / lfl is perpen¬ 
dicular to the orbit plane (opposite to the orbital 
angular momentum vector) I. 0 is the other corre¬ 
sponding co-ordinate of the right-handed system. A 
range angle 0 T is defined as the angle between and 
the radius vector at the estimated injection point on 
the orbit. 

Orbit 0rbi > 



meridian 
at launch 

Fig. 1. Orbit plane geocentric inertial co-ordinates. 


Closed-loop guidance 


121 


An estimation of the injection state on the orbit can 
be obtained using the orbital parameters and the 
estimate of the range angle d T as: 

V T = VolL sin(y D - 0 T ) + 4cos(y D - 0 T )], 

Rt = Rd [4 cos (Oj ) + 4 sin(0 T )], (1) 

where V D ,R D are magnitudes of velocity and radius 
vectors at specified true anomaly, and y D is the flight 
path angle corresponding to the desired true anomaly 
at injection. 

In order to determine the thrust attitude in real 
time, the following simplifying assumptions are 
made: 

(1) CLG begins after a short finite period of launch 
during which, the open loop guidance steers the 
vehicle beyond land mass constraints and dense 
atmosphere. 

(2) An equivalent constant/uniform gravity vector 
representation for the vehicle moving around spheri¬ 
cal earth is used. 

(3) For the remaining part of the flight, nominal 
parameters of the vehicle defined a priori , are used in 
guidance computations. 

The equations of motion for a launch vehicle can 
be expressed as: 

y = (F(t)/m(t))u(t) + g, r = v, (2) 

where v and r are the 3-dimensional velocity and 
position vectors, u(t) is a unit vector along the 
direction of thrust. F(t), m(t) are the instantaneous 
thrust magnitude and mass respectively, g is a con¬ 
stant acceleration due to gravity. Since u(t) is a unit 
vector it satisfies the constraint, 

«(/)• * 3 ( 0 = 1 - ( 3 ) 

The vector u(t), and the duration of thrust are 
decided at each instant from the information on 
current position and velocity vectors from INS (Iner¬ 
tial Navigation System), the predefined target and the 
fuel optimality. For a rocket with non-throttlable 
engines, fuel optimality implies minimum time of 
burn. The associated performance index can be given 
as, 

f' f 

J = dt. (4) 

Jo 

The necessary conditions include, state and co-state 
equations and control constraint relations [11,12] as 
follows: 

4 = 0, 4 = -P r (5) 

where P,, P t , are the co-state vectors corresponding to 
r and v respectively. On solving eqn (5) one gets, 

P, = constant, P t , = - P r t + C. (6) 

The transversality conditions involve certain state 
variables being specified at the terminal point. The 
remaining adjoint variables are set to zero at the 


terminal time, since, the performance index does not 
contain a penalty on the unspecified state variables. 
The control constraint relation leads to, 

(F(t)/m(t))P v = 2P u m, (?) 

where P u is the adjoint variable associated with the 
constraint in eqn (3). Therefore, the unit thrust vector 
u(t) lies in a plane defined by P, and C, [eqns (6-7)] 
even for the 3-dimensional launch trajectory as it is 
along 4,. 

In general, P,, P„ can be represented in any inertial 
frame of reference. Since, the thrust vector lies in a 
plane throughout its operation, an appropriate guid¬ 
ance frame is defined as shown in Fig. 2. The 
coordinate 4 is defined along — P r , 4 is perpendicu¬ 
lar to the 4 > n the ( 4 , 4) plane, and towards P,. P,, 
can now be resolved along I xc and I yc as: 

P v =C x [L + fycP(t-T c )], (8) 

where C = f xc C x + I yc C y , P = P,/C x and T c = C y / 
P,. Since, u(t) and P„ are collinear, one gets, 

U(I) = [4 + IycP(t - 4)]/[l + P\t - TV) 2 ] 1 ' 2 . (9) 

The optimal thrust vector lies in the 4> 4 pl ane - 
Therefore, the changes in position and velocity vec¬ 
tors due to thrusting alone necessarily lie in this plane 
(V thR , R thR ). Hence, this plane can be called as correc¬ 
tion plane. The incremental position vector can be 
split along V thR (i.e. R lhm ) and a direction perpendicu¬ 
lar to V thR (i.e. R thn ) in the correction plane (Fig. 2). 
For launching the payload into the orbit, there is a 
freedom in terms of injection on the orbit. This 
enables the R thm component to be free, since for every 
choice of R thm , there exists a possibility of achieving 
the orbit injection requirements. The other com¬ 
ponent R lb „ gets suitably redefined for each choice of 
R, hm . Now, the costate vector P r can be split along 
Rthm and R lhn as: 


4 = unitjR.hJ.P™ + unit(R thn )P rn . 


4 




ult). 

A 

-P, 



I»C 



X 

\c 


^thn 






N 

Fig. 2. Correction plane. 
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Since, R, hm is unspecified, the corresponding transver- 
sality condition is, 

P m (ti) = o. 

Since, P r is a constant vector, its component P rm is 
identically zero. Therefore, — P m is collinear with ! yc . 
and hence, V lhR which is along R thm is collinear with 
. The use of a correction plane helps to resolve the 
contribution due to thrust on position and velocity 
vectors. From Fig. 2, the optimal steering law is given 
by the linear tangent law of the form, 

tan W(t)] = P(t-T c ). (10) 

Equation (10) is an important milestone in the devel¬ 
opment of CLG. The guidance parameters P, T c , f xc 
and t yc in eqn (9) are determined at each instant of 
guidance computation from the knowledge of the 
present and target positions and velocities, the contri¬ 
butions due to thrust and gravity. After determining 
above mentioned parameters, the pitch and yaw 
commands ii j/ y (t) and their time rates can be 
computed with respect to LPI coordinates by project¬ 
ing u(t) on pitch and yaw planes as: 

<P P U) = tan '[(I xcx + I ycx P(t - T y ))l(I xc; + ff ; P(t - r £ ))j 

* y (0 = tan+ I vcy P(i - TJ)/K], 

<Ap(0 = P[(I X cJyP X - 1^)1 K 2 ], 

MO = + Iy*nt ~ T e )))l 

(1 + P\t - T c ) 2 )], (11) 

where, I KX , I KJ , I XC! , I ycx , /,„ are x,y, z components of 
unit vectors I K , I yc in LPI coordinate frame respect¬ 
ively, and 

K 2 = [(/„ + I ya P(t - T c )f + (I xcx + I ycx P(t - T c )) 2 ], 
k = (P/K)[P(t - r e )(l - i; c> .) - I X M- 


3. IMPLEMENTATION OF THE GUIDANCE LOGIC 


To find guidance parameters mentioned above, 
estimation of velocity and position increments to be 
gained through thrusting, are required. These re¬ 
quirements can be calculated after estimating the 
effects of gravity during the future trajectory, and 
finding a suitable injection point on the orbit. On the 
other hand, a proper choice of guidance parameters 
fixes the thrust profile. The corresponding velocity 
and position increments are then matched with the 
above requirements. 


3.1. Calculation of thrust effects 

The changes in velocity and position due to avail¬ 
able thrust force F(t) acting along optimal u(t) 
direction can be written as: 

'Tgo 

(F(i)/m(f))fi(f)dt, 


*3 

Jo 


PlhA — 


Tgo 


(F(s)/m(s))u(s)As 


dr, (12) 


where V IhA , R thA are the velocity and position vector 
increments though thrusting, and Tgo is the total 
duration of thrust. These integrals can be evaluated 
exactly as explained later in Section 3.1.1. One can 
also use a simpler, but approximate method of evalu¬ 
ation of these integrals using a series expansion for 
the denominator of u(t) in eqn (9) as: 

[1 +{P(t - T c )) 2 ]~' 12 = 1 - (\/2)(P(t - 

+ (3/8)(P(r - T c )) 4 --.... (13) 

This series converges for P(t — 7j.) < 1, which implies 
that eqn (13) is valid for i/f(r)<45°, even though, 
optimal steering law is valid for \j/(t) <90°, as seen 
from eqn (10). Using this series as a polynomial in 
time variable t, V lhA , R lhA can be expressed in terms 
of integrals of the form: 

[■Tgoj 

Isnj= (Fj(t)lirij(t))t" d/, 


fTgcyf r> 

Idrij — (F / (s)/m,(s))s" d.s] dl. (14) 

Jo LJo 

where, index j corresponds to a particular stage of the 
multistage vehicle. Tgo y is the remaining burn time of 
yth stage. ldn t can further be expressed in terms of 
single integrals, on integrating by parts as, 


Idrij = Tgoj 



(Fj(t)lmj(t))t n At 


*7*0/ 

0 


(F,{t)lm,{t))t'" 


df, 


= TgOj Isn, - Is(n + 1),. (15) 

Equations (14) and (15) are valid for any thrust 
model. Most of the liquid propellant rocket engines 
can be modelled as: 


(F l (t)/m J (t)) = [Fj(t)lrhj(t)]l(TaUj - t)), 
= VeXjUTau, — t). 


where m 0j = vehicle mass at present instant t„, 
riij = constant mass flow rate of j th stage, 
VeXj— Fj(t)/mj(t) is a constant exhaust velocity for 
j th stage, and Tauj = m 0j lriij is a vehicle time constant 
updated at each guidance cycle when m 0j is changing 
with t 0 . With this thrust model, the integrals in eqn 
(14) become [7, 11,12]: 


I.m l = 


' Tgoj 

[VeXjKTau,■- t)]t” At. 
o 


On further simplification, the above integrals can be 
expressed as. 


Isnj= -(1 /n)Tgo" Vexj+ Taujls(n - 1),, (16) 

with IsOj = VeXj\n[TaUjl{TaUj — t)]. Using eqn (16), 
all thrust integrals can be found sequentially starting 
with IsOj. These integrals are evaluated at every 
guidance cycle instant, only for the current stage 
(j = cs) and the final stage. The integrals for inter¬ 
mediate stages need to be computed only once at the 
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beginning of the closed-loop guidance, as these do 
not vary during cycle to cycle computations. 

Neglecting third order and higher terms in bino¬ 
mial expansion in eqn (13), relations (12) can be 
expressed for an N --stage vehicle as. 


V thAx f* 


0 = L 


RthA ~ f 


N 


I (VAX lJ+ VAX 2J ) 


L; = cs 
N 


Z(VAY v +VAY v ) =V thAy , 


J m cs 


'} 

']■’ 


£ (RAXy+RAX: 


2 /) 


+ 7„ 


L/-cs 


/) - 


where suffixes x and y indicate I x , 7 VC components, 
and, 


VAXij = 7s 0,, 

VAX y = -(7> 2 /2)[7s2y - 27s 1, T CJ + 7s0, T 2 ], 
My ij = 7sl J -/sO J 7' c „ 

FT y 2 , = - (7> 2 /2)[7s3y - 37s2y T CJ 
-2Is\jTl + IsQjTlj], 

RAX U = 7(70, + r r> FTAjy, 

= - (P 2 /2) [Id2j - 2Id 1,7’ c/ 

+ 7c/0 ; r 2 ] + t rj VAXy, 

RA Y v = Idlj+ t rJ VAXy, 

RA Y 2j = — (7* 2 /2)[7<73y - 37*72, 7; + 37*71 y T 2 
+ 7*70, T\j\ + t rJ VAX 2j . 


Here, T cj =T C — T sj , and T sj is the starting time of yth 
stage with respect to present instant t 0 , except for 
j = cs, for which it is zero. The remaining flight time 
till the injection, after the completion of yth stage, is 
represented by t rJ . 

3.1.1. Exact evaluation of thrust integrals. For eval¬ 
uating the thrust integrals in eqn (12) exactly, con¬ 
stant thrust, uniform mass flow rate model as 
explained earlier is used. The effects of thrust can be 
expressed in terms of the following integrals defined 
for the yth stage as: 


4 = 


VeXji[(T Ul - t)R\ 2 ] dr, 


4 = 


[VeXjIR'f] df, 


(17) 


where Tu, = Tau t + T sj , T tj = T sj + TgOj, R, = 1 + 
P\t — T c ) 2 . R, can be expressed as a polynomial in t 
as: 


R, = c,t 2 + b,t + a,, 

where c, = P\ b, = -2P 2 T C> a, = 1 + T\P 2 . On sub¬ 
stituting a new variable Xj defined as Xj = l/(Tuj — t), 


I xJ can be expressed in a similar manner to I tj in eqn 

(17) as: 


7 XJ = J (VeXj/R'J 2 ) dx,. 

J x jl 

where R xj = c xj t 2 + b xj t + a xj , with c xj = 1 + P 2 (Tu ] — 
T c \ b xj = —2P\T Uj — T c ), a xj =P\Tuj-T c ) and 
x n =\/(Tuj-T S j), x j2 =l/(Tuj-T eJ ). Evidently, c, 
and c xj are always positive, they being quadratic 
functions. Hence, the above standard definite inte¬ 
grals can always be expressed after performing inte¬ 
gration, as: 


Ixj— (VeXj/(c xj y 12 In 


2(c„ R y ) li2 + 2c,, X: + b 


XJ « •'XJ 

T', 


X J 2 


Iij—(.Vexjl(c,) 112 ) ln[2(c, R,) >l2 + 2c,t + b,]. (18) 


With the above definitions, relations (12) become. 


N 

v = i Y r 

T thAx *xc C* *xj 

j = cs 

0 = LP £ ((Tu xj -T c )l xj -I,j) 

j= cs 


R,hAx = LlTgo V thAx + £ (7,y - I XJ Tuj)\ 

J m cs 

RthAy 7 ;c [ TgO F t j, A y 

+ /> £ [VexfiR.^lc.) 

i- « 


-(Tu,-T c - (b,/2c,))I,j 
-Tuj(T Uj -T c )I xJ ]]. 


(19) 


3.2. Computation of gravity effects 

In order to obtain accurate injection, an equivalent 
uniform gravity field over a flat Earth has been 
defined [9-12]. Two separate gravity vectors g p and g v 
are defined such that change in position and velocity 
along the guided trajectory over a spherical Earth 
(inverse square field) is approximated to that over the 
flat Earth as follows, 



“ 

r Tgo 

ir r Tgo -i 

gv = 


W|r 3 (OI)r(/)dr f 

dr 


— J 

o J/ 

Jo 


Tgo 


(p/\r\s)\)T(s)ds dt 


1/ 

■ r Tgo 

* t 

ds df 

J/ 

-Jo . 

0 


( 20 ) 


The changes in velocity and position vectors (V g , R g ) 
due to gravity are: 

V g = Tgo g v , R g = (l/2)7go 2 g p . (21) 

The integral terms on the numerator of eqn (20) can 
not be evaluated directly as r (t) depends both on 
gravity and thrust acting simultaneously. Hence, V g , 
R g are found using Enke’s method [14], In this 
method an osculating reference trajectory is estab¬ 
lished, using only central force motion. The incre¬ 
ments in the state variables (position and velocity) 


AA 25/3—B 
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due to thrust are added to the propagated state on the 
osculating trajectory at the rectification points. Since, 
thrust integrals are known for individual stages, the 
burnout instants of each stage can be the natural 
rectification points. On propagating state from the 
time t 0 for a time step of A t = t l ~t 0 , using an 
osculating trajectory one gets r 1; at time t ] (Fig. 5). 
The correction <5 r is defined as the difference between 
actual position vector r at r, and propagated position 
vector r, , i.e. 

<5r(/,) = r(t,) — r,. 

The <5 r can be thought of as resulting from accelera¬ 
tions due to perturbing forces along the osculating 
trajectory. Such an acceleration can be given as, 

<5r(ti) = f(/i) —r,. 

At all rectification points, the initial states of the next 
osculating trajectory are corrected to the true state of 
the vehicle. Hence, <5r is zero at these instants. It 
increases gradually as one moves along the osculating 
trajectory. At each rectification point, the correction 
<5 r is computed in stages, i.e. <5r = 5r, + <5r 2 • ■ •. The 
first correction <5r,, which is dominant, is due to the 
application of the thrust vector for the duration At. 
In reality the gravitational and thrust forces act 
simultaneously. The perturbations <5r 2 , <5r 3 , • • • due 
to coupling between the above two forces are com¬ 
puted sequentially [14] until the subsequent correc¬ 
tions, <5r,, are small. The oblateness effects can also 
be considered, if very high accuracy is desirable. The 
state variables are similarly obtained at all rectifica¬ 
tion points until the thrust cutoff where the state 
estimates (i.e. V Te ,R Te ) thus obtained are generally 
quite close to the desired target conditions. Then, g v . 
g p , are computed as, 

gv = (V Te — V 0 — V thA )/7go. 

gp = 2(R Te - r 0 - V 0 Tgo - R thA )/7go 2 . (22) 

3.3. Solution to guidance parameters 

The solution to the guidance problem consists of 
finding the parameters that decide the thrust steering 
angles and the time of thrust cutoff such that follow¬ 
ing constraints are satisfied: 

^thR = V lhAx , 

0 = V thAy > 

RthR Acc Ax A v f^thAv ■ (23) 

where the required velocity and position V thR , R thR 
are defined as, 

V th R = V T — V 0 — g v Tgo, 

RthR = Rt - Ro - (i/2)g p Tgo 2 - Tgo V 0 , (24) 

The right-hand side of eqn (24) depends on values of 
P, T c , Tgo, i xc and I vc . The gravity vectors g v , g p , in 
eqn (24) also depend on all these parameters. V T , R T 


depend on Tgo and 0 T . Solving eqns (23) and (24) for 
the guidance parameters is very cumbersome, as these 
equations are highly nonlinear. These parameters do 
not vary by a large quantity from one guidance cycle 
to the next. Therefore, it is possible to estimate 
gravity contributions fairly accurately from an earlier 
knowledge of the parameters, since the sensitivity of 
such an estimate is not large. The guidance par¬ 
ameters are refined using differential corrections 
while holding g v , g p constant. However, at the in¬ 
itiation of the closed-loop guidance, one may not 
have reasonable estimate of parameters to compute 
g v , g p to the requisite accuracy. Therefore, the above 
two-step procedure is repeated. 

3.3.1. Partial differential corrector. The vectorial 
constraints (25) are reformulated as scalar relations 


for the sake of simplicity as. 


' ' ^AhR ^thR “ ^AhAx ’ ^AhAx 

(25a) 

Tl = RthR ‘ RthR ~ RthA ‘ R|hA 

(25b) 

T} = V thR • R thR — V thA • R, hA 

(25c) 

T4 ^thAy ^AhAy • 

(25d) 


By the proper choice of Tgo, P , and T c . F, to F„ can 
be driven to zero. T c is initially solved using the 
expression for F t . Then, Tgo, P, 0 T are solved simul¬ 
taneously by a differential corrector method using 
eqns (25a-c). With these, V, hR and R, hR are updated, 
and Ac, Ac • are computed using, 

Ac = Unit vector (V thR ), 

Ac = Unit vector (R thR x/jx Ac . (26) 

Finally, the corrections to T c as a result of changes in 
Tgo, and P, is provided by the relation, 

dT c = ( dTJSTgo ) ATgo + (dTJdP) AP. (27) 

Here, partial derivatives of T c are obtained from 
relation in eqn (25d). The required partial derivatives 
are evaluated only once in a guidance cycle except 
during the first cycle on ignition of each stage. This 
reduces guidance computations substantially without 
sacrificing accuracy or optimality. This is so because, 
subsequent corrections are very small to cause ap¬ 
preciable changes in partial derivatives. 

3.3.2. Starting estimates of guidance parameters. 
Partial differential corrector methods described above 
need good starting estimates of the guidance par¬ 
ameters for fast convergence. Only the Tgo estimate 
based on ground simulations needs to be supplied 
externally. All other parameters can be estimated 
onboard at the starting of the CLG. However, for 
each subsequent guidance cycles, simple starting up¬ 
date formulae can be derived using the parameters 
from the previous cycle. 

Initial estimates at the start of CLG: The first 
approximation to range angle d T is obtained by 
averaging the angular velocity of the vehicle at the 
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guidance initiation time T a and the orbital angular 
velocity at injection, using Tgo estimate, as, 

0 t = (1/2)[(V o -/ zO )/(R o - 4 o ) 

+(V D cos y D )l(R D ))Tgo. (28) 

Here, V Q , y D , R D can be computed from the orbital 
parameters. / rfl , can then be calculated with 
the additional knowledge of position and velocity 
vectors at T 0 . Now V T ,R T can be found using 
relations in eqn (1). Using the present position and 
estimated injection position, the gravity related vec¬ 
tors are computed from a two-point approximation 
as, 

g. = (l/2)( g(R 0 ) + g(R T )), 

g p = (l/3)(2g(R 0 ) + g(R T )). (29) 

From the above approximations, V thR , R thR are calcu¬ 
lated [eqn (24)] and in turn /„, f vc are found as 
explained earlier [eqn (26)]. Then considering only the 
first term in the binomial expansion (eqn. (13)), 
estimates of T c and P are obtained as, 

i (Is\j+T SJ IsOj)jl^ £ A0,j 
P = — (R th R ' / vc )/1 — IdOj(Ty + Is l ; //sO ; )) 

j - cs 

+ t^/s\j + IsOj(T^ + h\jlhOj))]. 

Starting approximations at guidance cycles: 

At subsequent guidance cycles, the first estimates 
of guidance parameters are made from their values in 
previous cycles. Tgo and 0 T at the present instant are 
given by, 

Tgo — Tgo {M) — At, 9j = — A0 T , 

where At is the guidance cycle interval, and A0 T is the 
angle between projections of present position and 
previous position on the orbit plane. The value of P 
is retained as in the previous cycle. Gravity related 
vectors are also obtained similarly, by, 

8v [Svtoid) Tgo (a]d) A ]/ Tgo, 

8 P = 2[(g P (oid> Tgo\ 0 | d ,)/2 - AV g Tgo - AR g ]/Tgo 2 , 
where 

AV g = (l/2)[g(/J 0(old) ) + g(R 0 )](A/) 

AR g = (l/6)[2g(/? 0(old) ) + g(/? 0 )](Ar 2 ). 

The guidance scheme is very sensitive to gravity 
estimates. Depending on injection accuracy required 
AV g , AR g can be found by orbital propagation back¬ 
ward in time for the guidance cycle interval. But this 
leads to increased computational load. 

3,3.3. Terminal computations. As the launch ve¬ 
hicle approaches injection point, very large turn rates 
are demanded if all the CLG computations are 
carried out until the end. This tendency towards 


instability results from the vehicle modelling errors, 
and the demand of the zero injection error to be met 
in a very short time to go (Tgo). In order to overcome 
such undesirable effects, guidance parameter updates 
are stopped well before thrust cut off and previously 
computed parameters are used till the end. Only the 
velocity magnitude constraint equation [eqn (25a)) is 
solved to find the Tgo accurately. The faster updates 
of Tgo using simpler computations are necessary as 
the injection accuracy in terms of eccentricity depends 
greatly on accurate thrust cutoff. 


4. SIMULATION STUDY 

The simulation study of the CLG logic is made, 
using typical launch vehicle data (Table 1). The 
launch vehicle consists of four stages with a relatively 
long coast phase before the ignition of the final stage. 
The position and velocity vectors at the start of CLG 
are also given in Table 1. The parameters of the 
desired 900 km sun-synchronous orbit are given in 
Table 2. 

The average constant thrust and constant mass 
flow rate model of each stage in the vehicle are used 
in CLG computations. Simulation of the dynamical 
system is made using Adam's Bashforth method, with 
a step interval of 0.1 s. The guidance scheme consists 
of major and minor cycles. The duration of each 
major cycle is 4 s and it consists of four minor cycles 
of 1 s each. At the beginning of each major cycle all 
guidance parameters are recomputed. The pitch and 
yaw angles and their rates [eqns (11)] at each minor 
cycle are computed using these parameters. In the 
terminal phase the guidance computations are sim¬ 
plified with the updating of only Tgo at a major cycle 
of 0.5 s. The terminal phase starts when Tgo < 5 s. 
The thrust integrals are computed using the series 


Table 1. Vehicle characteristics and its state at CLG initiation 


CLG starting time = 204.4 s 
(with respect to 
launch instant) 

Tgo estimate at = 795 s 

CLG start-up 

Vehicle state at Position (km) 

(LPI) CLG start-up 

X = 285.96 Y = 

Velocity (km/s) 

AT = 2.51 Y = 

12.96 

-0.89 

Z = 6550.4 

Z = 1.10 

2nd stage 

3rd stage 

4th stage 

Launcy vehicle data (liquid) 

(solid) 

(liquid) 

(a) Ignition time (s) -104.4 

50.1 

481.9 

(with respect to 



CLG initiation) 



(b) Total burn time (s) 153.5 

91.75 

400.0 

(c) Exhaust vel (m/s ! ) 2860.3 

2857.4 

3017.7 

(d) Initial mass (kg) 54589.6 

11908 

3696 

(e) Mass-flow rate (kg/s) 244.3 

79.5 

4.6 

coast 1 


coast 2 

(a) Initiation time (s) 49.1 


141.9 

(with respect to 



CLG initiation) 



(b) Duration (s) 1 


340 
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Table 2. Desired and achieved target conditions 


Target conditions Last stage 



(km) 

e 

(x 10~ 7 ) 

i 

(deg) 

n 

(deg) 

Altitude 

(km) 

Velocity 

(km/s) 

burn time 

(s) 

Specified 

7278.137 

0 

98.99 

263.76 

900 

7.400 


orbit 

Case (1) 

7278.140 

3.9312 

98.99 

263.76 

900.001 

7.40046 

383.78 

Case (2) 

Thrust Hi 

7278.173 

612.23 

98.99 

263.76 

900.000 

7.40048 

374.28 

Lo. 

7278.127 

574.74 

98.99 

263.76 

900.001 

7.40045 

393.59 

Case (3) 

7278.140 

2.5266 

98.99 

263.76 

900.001 

7.40046 

383.32 

Case (4) 

7278.140 

2.2451 

98.99 

263.76 

900.002 

7.40046 

383.67 

Case (5) 

(1) Series 

integrals 

7278.140 

2.1465 

98.99 

263.76 

900.002 

7.40046 

383.94 

(2) Exact 

integrals 

7278.140 

2.1465 

98.99 

263.76 

900.002 

7.40046 

383.94 



Fig. 3. Pitch and yaw vs time (nominal thrust). 


method in cases 1-3, and analytical formulations in 
case 4 of the simulation studies listed below: 

(1) Simulation model of the vehicle is identical to 
the realistic data. 

(2) Off-nominal simulations use uniformly 1% 
higher or lower thrust value in simulation model 
alone. 

(3) Computational decay effects are simulated. 
The delay is assumed to be one second uniformly till 
the terminal phase. At the terminal phase delay is 
assumed to be negligible. 

(4) Simulation data are similar to case (1). 

Simulation results indicate, that very high accuracy 
of injection is possible. The results in the form of 
demanded pitch and yaw angles, the injection accu¬ 
racy, are given in Figs 3-7, and Table 2, for all the 



Fig. 4. Pitch and yaw vs time (exact integrals). 



100 300 . 500 700 900 


Seconds 

Fig. 5. Pitch and yaw vs time (1% high thrust). 
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100 300 500 700 900 

Seconds 

Fig. 6. Pitch and yaw vs time (1% low thrust). 



Fig. 7. Pitch and yaw vs time (1 s delay). 
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cases of simulations done above. The altitude and 
velocity profiles and all the guidance parameters are 
also plotted (Figs 10-13). The thrust steering angles 
in pitch and yaw vary smoothly over the entire 
trajectory with the exception of a short period at the 
end of penultimate solid propellant stage. The sub¬ 
stantial difference between the average thrust model 
used in CLG computations and the actual solid 
propellant model used in simulation is responsible for 
such large variations in pitch and yaw. It can be seen 
from Figs 3 and 4, that such variations are small if 
analytical formulation of thrust integrals are used 
instead of series approximation method. This is so 
because, the errors in the series approximation is very 
high at these instants as i {/(t) will be close to 45° [eqn 
(10)] [when ip(t) is 45° the series expansion in eqn (13) 
diverges]. This difficulty can be overcome by using the 
guidance parameters from the previous cycles, when 
the demanded steering angles increase over the pre¬ 
vious commands, beyond a reasonable range. The 
simulation results using such a strategy are shown in 
Figs 8 and 9, and Table 2, as case (5). This conse¬ 
quently needs a small extra burn time. The better 
approach, however, is to use the exact thrust model 
in the CLG computations, at the cost of increased 
computational time. 

On-board computational time requirements are 
estimated by measuring the time required for the 
same CLG computations to run repeatedly a large 
number of times. The average computation time can 
be found for the guidance cycle computations at 



Fig. 8. Pitch and yaw vs time (frozen prmts, series 
integrals). 



Fig. 9. Pitch and yaw vs time (frozen prmts, exact 
integrals). 


Table 3. Computational requirements 


CLG type _Series integration_Analytical integration 


(1) Code memory (kbytes) 

(2) Data memory (kbytes) 

(3) Execution time 


31 

3.5 



32 

4 


Stage 

CLG cycle time (s) 

Execution time (ms) 
AT-286 AT-386 

No. of 
diff. 

corrector 

iterations 

Execution time (ms) 
AT-286 AT-386 

No. of 
diff. 

corrector 

iterations 


0.0 

(CLG start 
cycle—all 

996 

165.1 

5 + 4 

1682 

230 

6 + 4 

2nd 

computations 
done twice) 

8 

387 

64 

3 

518 

88 

4 


16 

361 

60 

3 

454 

77 

3 


20 

361 

60 

2 

389 

66 

2 


52 

(stage start 

361 

60 

3 

340 

58 

2 

3rd 

cycle) 

96 

309 

51 

2 

270 

45 

i 


482.82 
(stage start 

253 

42 

3 

302 

52 

3 

4th 

cycle) 

486.82 

208 

35 

2 

223 

38 

2 

494.82 

167 

28 

1 

189 

32 

1 


863.32 

18.3 

3.3 

1 

18.9 

3.3 

1 


(terminal 
phase cycle) 
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100 300 500 700 900 

Seconds 

Fig. 10. Altitude and velocity vs time (series integrals). 



100 300 500 700 900 

Seconds 

Fig. 11. Altitude and velocity vs time (exact integrals). 


various representative instants along the powered 
trajectory. The following machines are used to run 
the CLG computations. 

(1) IBM compatible (Prakruthi AT/386) PC/AT- 
386 with 80386 Microprocessor running at 33 MHz 
with a coprocessor 80387 running at 25 MHz. 

(2) IBM compatible (Siva AT) PC/AT-286 with 
80286 Microprocessor running at 12 MHz with 80287 
coprocessor running at 10 MHz. 


The guidance calculations are compiled in Turbo 
PascaB compiler. The typical compiled code and 
memory requirements are given in Table 3. Since the 
final implementation can use dedicated hardware and 
software, the computational requirements indicated 
here are very conservative. The results of these 
measurements are given in Table 3. 


5. CONCLUSIONS 

A closed-loop, explicit, highly accurate and near 
optimal guidance scheme is given. This scheme is 
capable of steering the launch vehicles for sophisti¬ 
cated missions, which require large pitch and yaw 
manoeuvres and long range of trajectories. The ex¬ 
plicit scheme requires estimation of effects due 
to gravity and thrust from the present time till 
injection along the guided trajectory for deter¬ 
mining the guidance parameters and subsequently the 
steering angles. The gravity effects are computed 
using Enke’s method. The thrust integrals can be 
evaluated using the series approximation as well as 
the analytical formulation. Since, the effects due 
to gravity and thrust are related to the guidance 
parameters, a sequential algorithm is developed 
to determine the guidance parameters. Here, the 
computational time required with the exact analytical 
formulation is larger by approx. 10—35% compared 
to that with the series approximation method 
that uses only the first two terms in the expansion. 
The accuracy of target injection on the desired 
orbit is high. The injection errors in the semimajor 
axis lies in the range of few meters and in ec¬ 
centricity of the order of 10 -7 —10“ 5 . The algorithm is 
computationally feasible for real time implemen¬ 
tation. The computational delay does not affect 
the performance much. The CLG computations 
can be performed in reasonable time duration 
(Table 3). 



100 200 300 550 650 750 850 950 
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Fig. 12. Guidance parameters vs time (series integrals). 
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Fig. 13. Guidance parameters vs time (exact integrals). 
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